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Abstract 

We present a complete theoretical description of atomic storage 
states in the multimode framework by including spatial coherence in 
atomic collective operators and atomic storage states. We show that 
atomic storage states are Dicke states with the maximum cooperation 
number. In some limits, a set of multimode atomic storage states 
has been established in correspondence with multimode Fock states of 
the electromagnetic field. This gives better understanding of both the 
quantum and coherent information of optical field can be preserved 
and recovered in ultracold medium. In this treatment, we discuss in 
detail both the adiabatic and dynamic transfer of quantum informa- 
tion between the field and the ultracold medium. 

PACS numbers: 42.50.Gy, 03.67.-a 



1 Introduction 

In the interaction between optical fields and atoms, a photon can be ab- 
sorbed by an atom and then the excited atom can re-emit a photon by either 
spontaneous or stimulated emission. In this process the atom stores the en- 
ergy of the field and releases it back to the field. Recently, theoretical and 
experimental studies have shown that both the quantum state and coherent 

* Mailing address 
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information of the field can be stored in an atomic medium A very 

recent experiment witnesses that a signal pulse of light can be stored in an 
ultracold collective of atoms for up to a millisecond [[J. 

The basic scheme for storage of light information is carried out by electro- 
magnetically induced transparency (EIT) ||10|,||11[1. N three-level atoms with 



one upper level \a) and two lower levels \b) and |c) interact resonantly with 
both the signal and the control beams. The weak signal beam and the strong 
control beam drive the atomic transitions \a) — \b) and \a) — \c), respectively. 
Early investigations have shown that EIT permits the propagation of the 
light signal through an otherwise opaque atomic medium and that the group 



velocity of the signal pulse is greatly reduced [^2[-p4|. In a recent experi- 
ment [|TJ, when the control beam is turned off, the signal pulse is stopped 
and stored in atomic medium. This effect can be understood in terms of 
dark states. Dark states are combination states of the photon state and the 
atomic storage state. The bosonic quasiparticles in the dark state are called 
polaritons. When the strength of the control field is changed adiabatically, 
both the quantum state and the coherent information transfer between the 
signal field and the collective atoms @,||]. 

Due to the fact that dark states are eigenstates of the EIT interaction 
Hamiltonian, the transfer of the quantum state between the field and the 
atomic ensemble can be quasi-stationary in the adiabatic limit. Furthermore, 
mapping and storage of quantum information of light in an atomic medium 
may occur in a dynamical process without forming dark states. For example, 
Ref. J7j has shown that in a general Raman interaction, with a large detuning 
to the intermediate level, the nonclassical features of the quantum field can 
be mapped onto the coherence of the lower atomic doublet, distributed over 
the atomic cloud. However, a very recent theoretical study has shown that, 
in the EIT model, an adiabatic change of the control field is not necessary, 
and even a fast switching of the control field can be used in the writing and 
reading quantum information of the signal field 

In all of these models, collective atoms play the role of quantum memory. 
The mechanism of storage of the quantum field in a medium is based on 
establishing atomic storage states (or atomic memory states) which record 
all the information of the field. It is not surprising that photons can be 
transferred to atomic excitations in transition interactions. Physically, the 
question is how collective atoms record the coherence of an electromagnetic 
field. In the early theoretical work |T3| , Dicke studied the coherence and 
cooperation effects in atomic ensemble. He defined collective atomic opera- 
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tors as the sum of all the individual atomic operators, retaining the angular 
momentum properties. The Dicke states are the eigenstates of the angular 
momentum operators. In Dicke's paper, he considered two cases: the gas 
volumes have dimensions either smaller or larger than the radiation wave- 
length. For the latter case, the spatial phase distribution of the field has 
been included into the collective atomic operators. 

In the pioneering work contributed by M. Fleischhauer and M. D. Lukin 
0)0)111)0! the theory of atomic storage states basically consists of a single- 
mode description. In this paper, we present a complete description of atomic 
storage states in a multimode framework. A preliminary version of the briefly 



analysis described in this paper is given in Ref. fL6[ . Similarly, as in the 
second case of Dicke's work, we incorporate spatial coherence into the collec- 
tive atomic lower and upper operators. We indicate that, in the multimode 
description, atomic collective operators can behave as multimode bosonic 
operators, under the conditions of low atomic excitation and appropriate ra- 
diation wavelength; which is much larger than the average interval of atoms 
and less than the propagation length of medium. Aside from the original 
definition of Dicke states, we introduce atomic storage states with explicit 
expressions by containing spatial coherence of the radiation field, and indi- 
cate that they are Dicke states with the maximum cooperation number. The 
significant advance is that, under the conditions shown above, a set of multi- 
mode atomic storage states is established in correspondence with multimode 
Fock states of the electromagnetic field. This gives a better understanding of 
how both quantum and coherent information, of electromagnetic fields, can 
be preserved in atomic media. A detailed theoretical description of multi- 
mode dark states in the EIT model is discussed. Furthermore, parallel to the 
"stationary polariton" in EIT, we show the "dynamic polariton" formed in 
coupled harmonic oscillators. This illustrates the mechanism for the dynamic 
quantum transfer between field and macroscopic matter. 



2 Atomic collective operators with the bosonic 
commutation 

We consider N ultracold collective atoms which are approximately stationary 
at their positions. At very low temperature close to the critical temperature 
for Bose-Einstein condensation!]], the average kinetic energy of atoms is 
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greatly reduced. On the other hand, at low temperature, atoms are densely 
packed within a wavelength of optical field. The free path of an atom is much 
less than the wavelength, hence it confines the range of atomic motion. What 
the "still atoms" means is that, in the characteristic time of the system, 
the scale of motion for the centre-of-mass of the atoms is much less than 
the wavelength of the optical electromagnetic field involved. The two levels 
of atoms | b) and |c) interact with some optical field of wavevector k. We 
assume that N is a large number and the largest proportion of the population 
of the atoms is in level \b) throughout the system evolution, so that the 
completeness relation is given by 

N N 

^ = E(I^I + I^I)-EI^>M C 1 ) 

3=1 3=1 

It is not necessary that \b) is the ground state, for instance, in the case of 
EIT the level |c) can be lower or equal to \b). For the sake of convenience, 
we call | b) the "ground" state and |c) the "excited" state. 

In the interaction between field and atomic medium, the spatial coherence 
of the field affects only the local atoms. In the approximation of " still atoms" , 
the j-th atom located at position Zj suffers a local field strength with a phase 
exp(ikzj) . For this reason, we define the lower and the upper operators of 
the collective atoms as 

1 N 

a k = —7=22\bj)(cj\exp(—ikzj), (2a) 
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N j=1 
1 N 



\cj) (bj\ exp(ikzj) , (2b) 



N j=i 



where k is the wavevector of the optical electromagnetic field interacting with 
the transition \b) — |c). We notice that this kind of collective atomic operator, 
containing spatial coherence, was first introduced by Dicke, who investigated 
the super-radiate effect in collective atoms, in the case of medium dimensions 



larger than radiation wavelength [15]. In order to avoid the difficulties of the 
occurrence of the center-of-mass motion of atoms which may destroy the 
coherence, in Dicke's paper, he assumed the molecules are so massive that 
their center-of-mass coordinates will be then treated as time-independent 
parameter in equation. Now, this assumption can be implemented in the 
development of the ultracold technique. 
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For the purpose of controllable storage, the atomic transition \b) — \c) is 
usually a multi-photon process including signal and control photons. Equa- 
tion pal) should be replaced by 



1 N 



a k = - 7 =^|6 i )(c i |exp[-i(A;-A; c )2; i ], (3a) 



1 N 



a k = —7=^2\bj)(cj\exp[—i(k + k c )zj], (3b) 



where k and k c are respectively the wavevectors for the signal field and the 
control field. Equation fl3a]) describes a Raman transition, whereas Eq. ( pT5| ) 
describes a two-photon cascade transition. 

The commutation relations for these atomic operators are written as 

[o-fc, o-fc'] = Ht, <?U = 0, (4a) 

N 

K,4] = (1/iV) ^(I^X^-I - |q,-><c,-|) exp[-<(fc - (4b) 

3=1 

The exact commutation of Eq. fl£5D for the same mode is readily obtained 

la k ,a{] = (l/N)J2(\b J )(b J \-\c,)(c J \). (5) 
i=i 

If N is a very large number and most of the atomic population rests in level |6) 
throughout evolution, by applying Eq. ([!]), then Eq. ( f£E| ) is approximately 
reduced to 

N 

[a k ,a{,] ~ (l/J\0£expH(* -*')*,]■ (6) 

Assuming that the atoms are in a string and the average interval of the 
adjacent atoms is d, which is much less than the optical wavelength, i.e. 
kd « 1, one obtains 

r-, i r-w ,\n 1 ~~ expUkNd] explikL] — 1 

X>xp[ !% ] = g«p|*b - 1)4 = 1 _ ex ; [M] l - iV LJ , ( 7 ) 

where L = Nd is the length of the atomic medium. However, this result 
is also true for a volume of atomic gas which is considered as a continuous 
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medium 

,, , f L N .., . , explikL] - 1 
exp[ikzj\ = j — exp[ikz\az = TV . (8) 

In the case that the length of the atomic medium is much larger than the 
optical wavelength, we obtain 

^1*1 = {o {k L>>l). < 9 > 

By applying the above result to Eq. (^), one obtains the bosonic commuta- 
tion relation for the collective atomic operators 

[er fc , o-J.,] ~ 6 hk >, ( 10 ) 

where we should assume (k—k')L >> 1, or, equivalently, A— A' >> A 2 /(2%L). 
For the parameters used in the experiment L — 339/xm and A = 589. 6nm, 
so that \ 2 /(2nL) w 0.163nm, Eq. (|T0|) is a good approximation for distin- 
guishable modes. 

We summarize the conditions for the collective atomic operators satisfying 
the multimode bosonic commutation as 

iV > n, (11a) 
A, L > d, (lib) 
AA/A > X/L, (11c) 

where n is the number of atomic excitations and AA is the mode interval. 
The low excitation limit ( p. lap has already been shown in the previous paper 
Q. The other two conditions imposed on the radiation wavelength, Eqs. 
( |llb| ) and ( |llc| ), assure the atomic ensemble containing and distinguishing 



the coherence, respectively. We will see in the next section that the atomic 
collective operators behave similarly to the creation and annihilation opera- 
tors of the electromagnetic field. 

3 Single-mode atomic storage states 

The "ground-level" state of the atoms can be compared with the "vacuum" 
state, symbolized in Ref. as 

\C°) EE l&x&a.-.&tf). (12) 
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When the single-mode creation operators of the collective atoms are applied 
to the "vacuum" state, one obtains 



N 



(4r|C°> = -i=(Ei c i)^l ex P(^)] \bib 2 ---b N ) (13) 



1 



E 1 c *i ' ' ' °in)(K ■ ■ ■ hjhh ■■■b N ) exp[ik(z h H h ^„)] 



E ■ • ■ c n • ■ • c *„ • ■ ■ exp[i/c(z il H h ^„)] 



where S{i n } designates that, in the summation, any two indices cannot be 
equal, because (\cj) (bj\) 2 \bj) = 0. We note that some states in the summation 
of Eq. fljjD , for which the sequence in the index set {i n } is exchanged, are 
the same and should be put together. For example, (i\ — 1, ii = 2, 13, ■ • • , i n ) 
and {i\ = 2,22 = 1, 23, • • • , i n ) represent the same state. For an ensemble 
{i n } of n elements, there are n! permutations which form the same state. By 
eliminating these repeated terms in the summation, Eq. ( |I3"D can be replaced 
by 



(4r|C°) = Y, "\h ■ ■ ■ Ci, ■ ■ ■ Q„ • • • b N ) exp[ik( Zil + ■■■ + z in )\, (14) 
where J2{i n } is defined as 



JV-n+lJV-n+2 N-l N 

E"= E E ••• E E- (15) 

{in} »l=l *2=2 i n -\=n-\ i n =n 

■* 

{i 1 <i2<---<i n ^i<i n } 

The summation of Eq. ([TJ) includes (J) = iV(iV - 1) • • • (JV - n + l)/n! 
terms. Now, we define a normalized atomic storage state 



^ " V N(N - 1) • U '.{N - n + 1) (16) 
x 2 • • • Cfc • • ■ q, n ■ ■ ■ b N ) exp[ik(z h H h z in )]. 

{in} 

Obviously, the atomic storage states with a different number of excitations 
are orthogonal to each other 

<c?|<T) = Km- (17) 
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In this definition, the superposition state of N collective atoms includes any 
possible combination of n atoms being in the level |c), while the corresponding 
spatial coherence is recorded in the phase of the wavefunction. Physically, 
it means that n photons can be stored by any combination of excited n 
atoms with an equal possibility, in correspondence with the nonlocality for 
photons. However, the coherent information of the field has been retained 
in the probability amplitudes. Note that the atomic storage state, with a 
definite wavevector k, is independent of position z, which disappears in the 
summation. 

By using definition ([16]), Eq. ([14]) becomes 



It is easy to check that 

= yi^v^+Ticr 1 ). (19) 

The above two equations are exact. However, in the limit N >> n, the 
corresponding approximate expressions are 



(4) n \c°) * y/riffl), (20) 

and 

a\\C n k ) ~ v^TTlCr 1 ). (21) 
The annihilation operator is applied to the "vacuum" state 

a k \C°) = 0. (22) 

In Appendix A, the general formula for the annihilation operator is proved 
as 



Equations (pJ|) and (^3|) give immediately 



a{a k \C n k ) = (l-^J:HC fc ")~n|Cn, (24a) 
a k a\\Cu) = (l-^)(n + l)|Cy>=i(n+l)|CJ). (24b) 



The approximations in Eqs. ( |2"3| ) and (|2~3| ) are valid in the limit N >> 
n. Equation (p4|) verifies again the bosonic commutation in this limit. If 
one admits both the bosonic commutation ( |I0D and Eq. fl2T|), by using the 
commutation 

[a k ,(alr]^n(a{r-\ (25) 

it can also obtain 

a k \C n k )^^\C n k ~ l ). (26) 
The atomic storage states are also the eigenstates of the population operators 

N 

ECfoXW) = {N-n)\CS), (27a) 

3=1 
N 

Ed^)(9lQ) = n\Ct). (27b) 

3=1 

(see Appendix A) 

According to Dicke's definition ( Eq. (47) in Ref. [15| ), the total angular 
momentum operators of the atomic ensemble can be described as 



Rkl 


= (VN/2)(*l + a h ), 


(28a) 


Rkl 


= (-zViV/2)(4-^), 


(28b) 


Rs 


= (N/2)(ala k -a k al), 


(28c) 


R 2 


= Rkl + Rk2 + -^3 


(28d) 




= (N/2)(a{a k + a k a[) + (N 2 /4)(a{a k - a k a\f. 





Using the exact relations of Eqs. (|T9|) and (p3|), one obtains 

R,\Cl) = i(2n-JV)|Cy>, (29a) 

R 2 \C n k ) = \N{\N + l)\C n k ). (29b) 

Therefore, the atomic storage state defined in Eq. (|I^) is the right Dicke 
state with the maximum cooperation number r = N/2. The discussion in 
this section exploits a new feature of Dicke state. The Dicke states with the 
maximum cooperation number play the role of number states in front of the 
collective lower and upper atomic operators. 
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4 Multimode atomic storage states 



The multimode case is concerned with how the information of the multi- 
mode photons is distributed in the local atomic excitations \cj). To see it, 
we firstly consider a simple case — the multimode single-excitation atomic 
storage state; that is, each mode contains only one excitation. We apply the 
multimode creation operators to the "vacuum" state 



L---4JC ) (30) 

I N N 

Y \ c h)(K\ exptykiZiJ ■ • ■ Y I c in ) (K I exp(ik n z in ) \ hb 2 ■■ -b N ) 



JV h=l i n =l 



Y '\h ■ ■ -c ix ■ ■ -c in ■ ■ -fc^exp^Mii H h k n z in )]. 



This equation is apparently different to Eq. ([U|) by the phase factors. 
Indeed, the exchanges of the indices in the summation contribute to the 
same atomic state, but, with different phase distributions. For example, 
{ii — 1, z 2 = 2, z 3 , • • • , i n ) and (ii — 2, z 2 = 1, 13, • ■ • , i n ) display the same state 
\cic 2 b 3 ■ • • Cj 3 • • • c in • • • 6jv) j but with the phase factors exp[i(kiZi + k 2 z 2 + 

hzi 3 H h k n z in )] and exp[i(kiz 2 + + k 3 z is H h k n z in )\, respectively. 

Mathematically, for a given atomic collective state \bi ■ ■ ■ q 1 • • ■ Q n • ■ • 6^), it 
allocates n\ phase factors due to n\ permutations for n elements. This means 
that an atom in the level \cj), located at position Zj, records the information 
of all the modes. Because the field is global, each atom in the medium expe- 
riences the field coherence of all the modes, and, vice versa, the field of each 
mode affects all the excited atoms. 
To simplify the sign, we define 

i k n} ■ {z in }i = {kiz h H h k n z in )i, (31) 

where {zt n }i stands for the l-th sequence of all the n\ permutations for n 
elements. Accordingly, Eq. ( |30"D can be written as 

7.! 



ViV {ln} , =1 

(32) 

where Y^'L n \ has already been defined in Eq. ([15]). In comparison with the 
single mode case, shown in Eq. (14), n-excitations in Eq. (|32D share the 
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phases of n modes. We define a multimode single-excitation atomic storage 

\Cl ■■■C l k )= 1 (33) 

" y/N(N-l)---(N-n + l) 

x "\h ■■■c il ---c in - ■ -b N ) ^exp[i{£; n } • {zjjj], 
{i n } 1=1 

which has been normalized, as shown in Appendix B. With the combination 
of Eqs. fl32|) and (|33|), we obtain 



< ■ -HO = p N - 1] - N !r- n + 1) n ■ ■ ■ clj - K ■ ■ ■ cij. 

(34) 

Now, we discuss the general case of multimode atomic storage states, s 
modes containing total n excitations can be generated by 



(4 1 ) mi ---(4j ms |C°> (35) 



where m\ + • ■ ■ + m s = n. By defining an index set as 

{^n } = 1 j * * * ) 'mi j ^mi+1 1 ' ' ' i ^rai i ' ' ' i ^n— m s +l j * ' ' j ^n) j (36) 

Eq. fl35|) can be written as 

(4 1 ) mi ---(4j m lC°) (37) 

^ / |6i---Q 1 ---Q n ---6 i v) 



x expfifciC^ + ■ • • + z imi ) + • • ■ + ik s {z in _ ms+1 + ■■■ + z in )\. 

Similarly, as in the previous cases, any particular atomic collective state 
|&i • • • ■ ■ • Cj n ■ ■ ■ b N ) is related to n\ terms in the summation throughout 
all indices. But, among these n! terms, the phase factor of each term will 
repeatedly appear mi! • • -m s ! times because exchanges of indices within a 
mode cause no difference. As a result, the remaining non- repeated phase 
factors terms are n\/{mi\ ■ ■ - raj). We define again 

{k ( r a) } ■ {ziji = (fctei + ■ ■ • + z i mi ) + ■■■ + k s {z in _ ms+1 + ■■■ + z in ) 



(38) 
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as one of these combinations with index I. Equation ( |T7D is written as 

(<7l) mi ---(4 a r°\C°) (39) 

I i n!/(mi!-m,!) 

= '"^^'-y." '>:■■■<■,. ■■■<■„■ ..b x} £ expt^M^R 
viv {in} l=1 

A general multimode atomic storage state can be defined as 

\Cg~.C£) (40) 
/ iri]! • • • m s \ 



V N(N -1)---(N -n + 1) 

n\ I (mi!'"OT 8 !) 

x £ > • • • ct, ■ ■ ■ Ci n ■ ■ ■ b N ) £ exp[z{^)} • {z in }l 

{in} 1 = 1 

which has been normalized (see Appendix B). Finally, Eq. (|39"D becomes 

(<r •••«)"* i (4i) 

* y/m 1 \---m a \\CZ 1 •■•(%'), 

where the approximation is valid in the limit iV >> n. In this limit, the 
creation of an excitation of mode I for a multimode storage state is written 
as 

< I^T • • • C? • • • C%) ~ v^T+ll^ 1 • • • C,7 +1 • • • C™;). (42) 
Similarly, as in the single mode case, by considering the bosonic commutation 



l0| ) and Eq. ([|l]), the annihilation of an excitation of mode / is written as 

**|CT • • • c% ■ ■ ■ cp) * v^iicz 1 ■ ■ ■ cTr 1 ■■■CT;)- (43) 



Equations (0), flU) and (gl|) imply that the state |C°) defined in Eq. 
([12] ) represents a vacuum state not only for a single mode but also for a 
multimode. Physically, it shows that the ultracold collective atoms are able 
to store a multimode field. As a result, we may derive 



"IK 1 ) * \C k2 C^), (44) 
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where (C^ 1 ) can be understood as either a single mode state or a multimode 
state with the vacuum for those modes other than mode k\. 

In the above theoretical description, we have shown that the atomic stor- 
age states are a duplicate of the Fock states of the electromagnetic field. The 



explicit expressions of atomic storage states Eqs. (jlq ) and fl40D display dual- 
ity of particle and coherence. The excitations may appear everywhere with 
an equal probability in the medium, corresponding with the nonlocality for 
photons. However, each excitation takes the local phase factors of the single- 
mode, or multimode, fields as the quantum probability amplitude, therefore 
recording the coherence. Consequently, it may establish the correspondence 
of two quantum system, the field and the atomic ensemble. This provides 
the basis for complete storage of quantum information of the bosonic field in 
an atomic medium. 



5 Dark states in EIT 

In the EIT configuration, the weak signal field interacting resonantly with 
the atomic transition \a) — \b) is described by the field operator 

E s (z,t) = (l/2)£ a(z)exp[i(k s z - uj s t)] + h.c. (45) 

= (l/2)£o X! a (l) exp[iqz] exp[i(k s z — u s t)] + h.c, 
q 

where £ is the field amplitude per photon and u s = ck s . c is the speed of 
light in a vacuum. The strong control field, resonantly driving the atomic 
transition \a) — |c) is assumed as classical 

E c (z, t) = — ^- exp[i(k c z - u c t)) + c.c, (46) 

where p ac is the dipole moment of the transition \a) — \c) and uj c = ck c . In 
the interaction picture, the interaction Hamiltonian is described as 

Hj = hJ2^ q a\q)a(q) (47) 

h N f 1 

~o H \ 9^2a{q)\aj){bj\ exp[i(k s + q)zj] + ft\aj)(cj\ exp[ik c Zj] + h.c. \ , 

1 j=l { q J 

where u q = cq is the detuning of mode q with respect to the resonant fre- 
quency uj s of the signal field. By defining the atomic collective operators 
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1 N 

Pab((l) = T7Hl a i)( & il eX P^(^ + ?) 2; i]' ( 48a ) 
iV 3=1 

1 N 

PaM = ttYI \aj)(c j \exp[i(k c + q)z j ], (48b) 
i=i 

the Hamiltonian fl4"T|) is written as 

Hj = h^ q a\q)a{q) - ^ a(g)p a6 (g) + fiiVp ac (0) + h.c.j . (49) 

The new quantum field operator defined in Ref. || is written as 

ip = cos 9a q — sin 9a q , (50) 
where 



cos# = n/y/Q? + g 2 N, sin6 = gVN/^ftt 2 + g 2 N. (51) 

The transition \b) — \c) concerns both the absorption of a signal photon and 
the emission of a driving photon. Replacing k by {k s + q) — k c in Eq. (0), 
one obtains the annihilation operator a q of the collective atoms in EIT 

N 

\bj)(cj\exp[-i(k s + q- k c )zj}. (52) 



1 



a 



Correspondingly, k should also be replaced by k s + q — k c in the atomic storage 
states. ip q satisfies the bosonian commutation relation as long as a q does 

bP q , 4>l>] = cos 2 9[a q , a\,} + sin 2 6{a q , a\] ~ 5 qq >. (53) 

It has been shown in Eq. (|5lD , that the parameter 9 is related to the strength 
of the control field. In the strong and weak limits of the control field, ip q tends 
to a q and a q , respectively. 

According to Ref. 0, the dark state is defined as 

\D:) = ^=(^ q r\o)\c°), (54) 

where |0) is the vacuum state of the signal field. The lowest dark state is 
designated as \D°) = |0)|C°). The quasi-particle in the dark state is called 
a polariton || . 
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Using Eq. ([H]), one obtains the exact expression of the dark state 



" n(n-l)---(n-m + l) N(N -!)■ ■ -{N - m+ 1) 

w - i^y) v ^ v 

xcos n - m #sin m %-m)|C7) (55) 



E 



ln(n - 1) • • • (n - m + 1) /iV(iV - 1) • • • (iV - m + 1) 



777 1 V /V r ' 

m=0 ' V ly 



X {W + ff>N) n/a q 

The dark states described above are orthogonal to each other since they have 
different quasiparticle numbers, but are not normalized. Under the condition 
N >> n, the dark state can be approximately written as 



\D n q ) ~ f(-l) w p ^ m" m + 1) cos "" m ^ sinm g|n-m)|C g m ) 



m=0 



The above expression of the dark state satisfies the normalized orthogonal 
relation 

(DS\D?) = 5 nm . (57) 

Eq. QSEp shows that when the parameter 9 is taken to be and tt/2, the 
summation in the dark states reduces to only the first and the last term 

\Dg) = \n)\C°) for 6 = 0, (58a) 
\D n q ) = (-l) n |0)|^) for 6 = n/2, (58b) 

respectively. Therefore, by varying 9 adiabatically, the quasi-particles can be 
transferred between the photon state and the atomic storage state. 
According to definition fl5"4|), one can obtain the exact expression 

rtw) = -^(i>ir +i \o)\c°) = v^+ri^r 1 ). (59) 

v n\ 

It is easy to check 

vg£>°) = o. (60) 
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As the same for the operator a q , with the help of the bosonic commutation 
relation (153T), one obtains for the dark state (|5T 



1> q \D») ~ MDT 1 )- (61) 
Moreover, the multimode dark state can be generated by 

K ■ -- D i s ) = v/ni , 1 .. nJ «r • • • (vino>ic°>. (62) 

They can be treated just like the multimode photon number states. 

In Appendix C, we have proven that, at the exact resonance, both the 
exact expression ([55]) and the approximate expression (|56| ) of the dark states 
are the eigenstates of the interaction Hamiltonian fl49[ ) with a null eigenvalue. 
A pulse of monochromatic light has a narrow bandwidth, and the detuning 
uj q from the carrier frequency u s is small. If we omit the first term in the 
Hamiltonian (|49|), the multimode dark states consisting of the pulse are the 
eigenstates of the interaction Hamiltonian. 

Assume that, at the initial time, a signal pulse is at a multimode state 

J2 a (Qir--,Qs)\ni---n s ), (63) 

while the cold collective atoms are, approximately, in the ground state \C°). 
The combined system of the signal field and the atoms is in the state 

|*(0)> = £ a( qi , ■■■,q s )\n 1 --- n s )\C°) = £ a(q u • • • , q a )\D% ■ ■ ■ D%) e=0 . 

{is} {q s } 

(64) 

When the control field is strong enough, the signal pulse can maintain and 
transmit through the medium. Note that |^(0)) is also the eigenstate of the 
interaction Hamiltonian with a null eigenvalue. If the control field is changed 
adiabatically to a very small level at a later time ti, the state of the system 
is also changed adiabatically to 

|¥(t0) = J2a( qi ,---,q s )\D^...D^) 0=n/2 (65) 

= E (-i) B1+ ~- H Mft, • • • , Q.M\C£ ■ ■ ■ c n q ;). 

{is} 

It forms an associate state for |\l/(0)). The whole of the quantum information 
of the signal pulse has been stored in the atomic medium, in the form of a 
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"negative copy", in which each excitation changes a 7r-phase. As soon as 
the control field returns to the previous level, the state fl6"4|) is recovered. 
Conversely, if Eq. (|65| ) is an initial state generated in other model, by turning 
on the control field, it will be converted to the corresponding optical field, 
enabling it to be seen. 

6 Dynamic quantum transfer in macroscopic 
matter 

Due to the fact that the dark states are eigenstates of the EIT interaction, 
quantum transfer processes between field and matter are quasi-stationary by 
adiabatically changing the control field. On the other hand, the transfer can 
be performed in a dynamic way, which has been described in the literature 
@, |§. In this section, we study a general description for dynamic transfer 
of quantum state between field and ultracold matter. The interaction con- 
figuration can be designed as, either, the parametric process or the Raman 
transition in which both a weak signal beam and a strong control beam 
interact resonantly with two levels of atoms. 

For simplicity, we consider a single-mode interaction. In the interaction 
picture, the effective interaction Hamiltonian is written as 

Hi = hn(aa j! + aV), (66) 

where fi is the Rabi frequency of the control beam, assumed as classical. The 
collective atomic operator o is defined by Eq. (^) and behaves boson-like in 
the low excitation limit. The model is well known as a coupled harmonic 
oscillator, and can be solved exactly. Here we illustrate this model again 
from a new viewpoint by introducing a very simple method for the exact 
solution of the state- vector evolution. For this model, it is easy to obtain the 
evolution of the operators in the Heisenberg picture, 

a(t)\ _ ( cos fit —i sin fit \ fa(0)\ 
a{t)) ~ \ -i sin fit cos fit ) \a(0))' { ' 

With this method, if the initial state can be written as |\I/(0)) = /(a(0), a(O))|0(O)), 
while the evolution of the state |0(O)) is already known to be |0(t)), we ob- 
tain 

= J7(t)|4r(0)) = ^(t)/(o(0),(r(0))|e(0)> (68) 
= U(t)f(a(V),o(Q))U-\t)U(t)\<d(Q)) = f( a (-t),a(-t))\Q(t)), 
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where U(t) = exp(—iHit/h). An initial Fock state for the signal pho- 
tons m and the atomic collective excitations n is represented by \m, C n ) = 
(1/ Vmlnl)^ (Q)] m [a^ (0)] n \0, C°). Since the evolution of the vacuum state, as 
|G(0)> in Eq. (|1), is known to be |9(t)> = exp(-i#/t/7i)|0, C°) = |0,C°), 
we obtain the evolution for an initial Fock state \m, C n ) 

\* mn (t)) = (l/VM)[a\-t)] m [a\-t)] n \0,C°), (69) 

where the dynamical operators are written as 

a f (-t) = C/'(t)o + C/'~ 1 (*) = a f cos fit - wx 1 " sin fit, (70a) 
<r*(-t) = U^U- 1 ^) = a f cos fit- io + sin fit, (70b) 

where = a^(0) and cr* = ^(O). According to Eq. (|69|) , the exact evolution 
for an arbitrary initial state S£mJ m 5C n ) * s therefore J2£ mn \^mn(t))- 

The state \^ mn (t)) described by Eq. ( |69D evolves periodically with the 
fundamental frequency fi and conserves the total particle number. For exam- 
ple, in the simplest case of only a single excitation |1, C°), the time evolution 
is |^io(t)) = cosfit|l,C°) — i sinfit|0, C 1 ). The entanglement between two 
subsystems increases in time and reaches a maximum at fit = 7r/4. Then, 
the entanglement decreases and the excitation transfers completely from one 
subsystem to another at fit = tt/2. In general, if the initial state has m signal 
photons and no excitation for atoms \m,C°), the time evolution is written 
as 

|*mo(*)) = (l/Vrn\)[ai(-t)] m \Q,C°) = (l/v 7 ^!)^ cosfit - ^ sinfit] m |0, C°) 

p Q v r.(m-j)\ 

The evolution states can be called "dynamic polaritons", since they are 
compatible with the "stationary polaritons" defined by Eq. (|56|). At the 
times fit = (l/2)7r, it and (3/2)n, the evolution state has been de-entangled 
to (-i) m |0,C m ), (-l) m |m,C°), and i m |0,C m ), respectively. This means, 
at certain times, the entanglement formed in the dynamical process can 
be cancelled and the excitations are transferred completely from one sub- 
system to another. When an initial state is an arbitrary superposition of 
Fock states for the signal beam and the "vacuum" for the ultracold atoms, 
i. e. |^(0)) = |$,C°) where |$) = J2 a m\ m ), h evolves to the states 
|0, \&-\C°), and |0,$W) at the times fit = (l/2)vr, vr and (3/2)tt, 
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respectively, and comes back to the original at Qt = 2ir. Here we define three 
states associated with |$) = ^a m \m) 

l$(±0) = Y,(±i) m <xmH, l$ ( ~ } > = E(- 1 ) m «-l^)- (72) 

Obviously, these associate states have the same particle distribution, but with 
different phase shifts in amplitudes. As a matter of fact, the phase factor 
in the superposition can be observed macroscopically. For instance, if |<&) is 
a coherent state \a), one obtains the coherent states again for the associate 
states with a particular phase shift, i. e. |$( ±J )) = | ± ia) and l*^ - - 1 ) = 
| — a). As for an arbitrary state, it is also true by means of the expectation 
value of the amplitude operator, such that ($ < - ± ^|a|$'- ± *- ) ) = ±i($|a|<3>) and 
($( - )|a|$ < - - ' ) ) = — ($|a|$). Though, in general, these associate states are not 
identical to the original, owing to a phase shift; the quantum information of 
the original state can still be faithfully preserved. It looks like a photograph 
and the corresponding negative copy. The associate states |<3>( ±l )) and |$ ( - _ - ) ) 
can be seen as "orthogonal" and "negative" copies of a quantum "picture" 
|$). 

In the dynamic quantum transfer, a quantum state of the signal field can 
be stored in and then retrieved from a medium by turning off and on the 
control field at a certain time. Because the dynamic polariton state defined 
by Eq. fl7I]) is the eigenstate of the free Hamiltonian, it will be preserved 
while the interaction is turning off. Similarly, the model can be extended 
to the multimode case provided the conditions for the radiation wavelength 
Eq. ( |TTD is satisfied. Thus, coherent information of the optical field can be 
transferred simultaneously. 

Finally, we indicate that, for a proper transfer in this model, collective 
atoms must be prepared in an atomic storage state. For ultracold matter, its 
initial state can be considered as the vacuum state \C°) approximately. If, at 
the initial time, the field is at an arbitrary state while collective atoms 
have been prepared in a superposition state |<3> 2 ) = P m \C m ), l^i, $2) will 
evolve to the corresponding states |$2 , $1 >$2 ) and |$ 2 j) ,$^) 

at the certain times mentioned above. It displays a complete swapping of 
quantum states for the coupled harmonic oscillators. 



19 



7 Conclusion 



In conclusion, we define collective atomic operators and atomic storage states 
by containing spatial coherence and illustrate the conditions under which the 
multimode collective atomic lower and upper operators are boson-like. We 
indicate the fact that the atomic storage states shown by definition ([II]) are 
Dicke states with the maximum cooperation number. The new feature for 
these Dicke states is that, in the low excitation limit for a large number 
of atoms, they behave as the Fock states of an electromagnetic field. The 
complete description and the deductive explicit expressions for the atomic 
storage states present better physical understanding of why the atomic en- 
semble can record fully the quantum information, both the excitation and 
the coherence, of an optical electromagnetic field. In addition to adiabatic 
quantum transfer by means of dark states in EIT, we discuss the mecha- 
nism of dynamic quantum transfer via dynamic polaritons which is formed 
by the fundamental interaction between the field and ultracold matter. A 
combination of adiabatic and dynamic schemes may find more applications 
in quantum information technology. 
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Appendix A 

First, we derive the exact equation ( p3]) . By applying definitions (0) and 
([T|), h gives 

°>\d) = ^|lMNexpHM ^... (jv "L„ +1) (AD 
x ^"|6 1 ---c il ---Ci n ---6 JV )exp[iA;(2; ll H V z in )\. 

{in} 

When the operator J^fLi \bi){ c i\ applies to a particular state |&i c$ n ••• 

it produces a superposition of n states in which |cj fe ) is sequentially replaced 
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by \bi k ). After this operation, Eq. (|7^) is a summation of n x rM = 
N(N — 1) • ■ • (N — n + l)/(n — 1)! states, in which n — 1 atoms are populated 
in the level |c). It is in fact N — n + 1 times \C^ 1 ). For example, a par- 
ticular state in |C^ _1 ), say, \c\ - ■ ■ c n _i& n ■ • • b^), comes from N—(n—l) states, 
|ci ■ • • c n _ic n 6 n+ i ■ • ■ &zv/j l c i ■ ■ ■ c n -ib n c n+ i ■ ■ ■ 6tv) ; •••> |ci • ■ • c n _i6 n 6 n+ i ■ • • Cat), 
in |C£). Thus, Eq. ( |73|) is written as 



= ;W ^-V-n + l) (Af -" + 1) <A2) 



x 51 • • • Cii • • • c i„-i • • • M exp[iA;(^ 1 + • • • + z^J] 

{in-l} 



(JV-n + 1) 



iV---(iV-n + 2), . 



iVV N---(N-n + ly J \ (n-l)\ 

V^ICT 1 ). 



■icr 1 ) 



'iV-n + l -.^ 



V N 

Then, we prove Eqs. (f27|). Equation fl27a|) is written as 

N 

E(l 6 <><«) (A3) 
i=i 



ni 



\ N---(N -n + 1) 

N 

x "\bi){bi\bi •••c il ---c in - • ■b N )exp[ik(z il H h O] 

Z=l {i n } 

N 

^2^"\bf ■ ■ c h - ■ ■ c in - ■ - b N ) 



x exp[ik(z h + ■■■ + z in )](l - 5 lh ) • • • (1 - 5 Un ). 

Because all the indices of ik are not equal to one another, one has 

(l-S Hl ) ■ ■ ■ (1-S lin ) = l-(S Ul +- ■ ■+S lin )+(5 lil 5 li2 +- ■ •)- ■ ■ = l-(S lh +- ■ -+5 liv 

(A4) 

Substituting Eq. (|73|) into Eq. (|73|) , one obtains 

N 

E(\ b i)( b i\ C k) = (N-n)\C n k ). (A5) 
i=i 
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Eq. (B7H) is proved as 



N 



E(l c ^» (A6) 
1=1 



nl 



V N---(N -n+1) 

N 

x E E "|c,)(ci|6i • • • • • • Ci n ■ ■ ■ b N ) exp[ik(z h H h 2jJ] 

z=i {i n > 



y N---(N -n + l) 

N 

x E E "\bi---c h ■■■c in - ■ -b N )exp[ik{z h + ■ ■ ■ + z in )](8 Ul + ■ ■ ■ + 8 Un ) 



1=1 {i n } 
n\C n k ). 



Appendix B 

In this appendix, we calculate the normalized coefficient of the multimode 
storage states. A multimode single-excitation state is defined as 

K ■■■ C ln) =aJ2"\bi---Cn--- c in •••ME exp[i{AvJ • {z in U (Bl) 

{in} 1 = 1 

The probability of each state in the above superposition is 



Eexp[i{/c n } • {ziji} 



1=1 



2 



^Mi{kn} ■ ({z in }i - {z in }j)] (B2) 
1,3 

nl 

n\ + Eexp[i{A; n } • {{z in }i - {z in }j)}. 



Then, we sum all these probabilities. For the first term of Eq. (|73|), it is 
simply 

„ nj = B|Jy(Jy -i)-( Jy -„ + i) = _ _ _ _ n + 

r 1 n! 

{in} 
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In the multimode case, one must find a mode with kj ^ 0. By using Eq. (Q), 
the summation J2 'L n \ to the second term of Eq. (|73| ) vanishes. Therefore, 
we obtain 

a = . (B4) 



the 



N(N -1)---(N -n + 1) 

Similarly, for a general multimode storage state defined by Eq. 
probability of finding a single state is 



n! / (mi!-m s !) 

£ exp[i{A«} " {ziM 
1=1 

n!/(mi!-m s !) 

E exp[i{A^)} • (K>, - 

n!/(mi!'-m s !) 



(B5) 



n! 



• • • mj 



+ £ exp[z{e s) }-({^}/-{^},)]- 



The summation to the first term of the above equation gives 



E 

{in} 



nl 



nii\ ■ ■ ■ m s \ 



(B6) 



n\ N(N -1)---(N -n+1) N(N - 1) • • • (JV - n + 1) 



mi! ■ ■ ■ m s \ 



nl 



mi! ■ ■ -m s \ 



With the same reason, the summation to the second term vanishes. The 
normalized coefficient is therefore 



a 



mil ■ ■ ■ m s \ 



N(N - 1) • ■ ■ (JV - n + 1) 



(B7) 



Appendix C 

We define a new collective atomic state in which n atoms are in the level |c) 
whereas one atom is in the level \a) 



K c n q ) 



nl 



^-i)-(i^ 1 £g" |6l - a, - Ci '- c -- M 



x exp[i(k s + q - k c )(z h H h z in )\ exp[i(k s + q)z{\ 



(CI) 
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where J2i^{i n ] designates the summation for index I which cannot be taken 
as ii, ■ ■ - i n . We have already indicated that, the state \C™) is a superposition 
of N(N — 1) ■ • • (N — n + l)/n! possible states \bi ■ ■ ■ • ■ • Q n • • ■ b^) in which 
n atoms are in the level |c) whereas the remaining N — n atoms are in 
the level \b). For one of these states, each of the N — n atoms being in 
the level \b) can be excited to the level \a). So the state \A q ,C q ) includes 
N(N — 1) ■ ■ ■ (N—n)/n\ such possible states |&i • • ■ a r ■ ■ • ■ ■ c in ■ ■ ■ b N ) with 
equal possibility. The state \A l ql C q ) has been normalized. The phase factor 
related to the excited atom I being in level \a) is exp[i(k s + q)zi], because the 
transition of the level \a) to the ground level \b) is connected with the signal 
field of the wavevector k s + q. State (|T3| ) can be obtained by the following 
operation 

N Pa M\c n q ) 

= (j2\ a i)( c i\ e Mik c zi]) \C n q ) (C2) 



\i=i / 

N 



V iV--(iV-n + l) g g >l) {Cllbl ' ' ' Ch ' ' ' ° in ' ' ' M 

x exp[z(/c s + g - H h O] exp[iA; c z/] 

V iV ft"h Ij ¥{in _ l}{in _ l} 

x exp[z(fc s + q - k c )(z, h H h Zi„_J] exp[i(fc s + g)z z ] 



ft! 



Ar... (i v_ n + i) 



V N---(N -n + l)\ (ft-1)! 

= v/ft^c;- 1 ), 

where the atomic operator p ac (0) has been defined in Eq. ( |48b|) . Similarly, 
we have 

N P M)\c n q ) 

= ^\a l )(b l \ex P [i(k s + q)z l ]j\C^) (C3) 

N 

E E Wi • • • cji • • • Ci n • • • b N ) 



ft! 



iV.-.(iV-ft + l)^ { ^ } 
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x exp[i(k s + q- k c )(z h H h z in )\ exp[i(k s + q)zi\ 



TV. 



Viv...(iv-nTT) / J w E / ,|6l --- ar -- Ql --- Ci "--" M 

x exp[i(k s + q- k c )(z h H h z in )) exp[i(k s + q)z t ] 



IV. 



V JV---(JV-n + l) 



,N - {N - n M^ 



nl 



The two interactions induced by two fields in the interaction Hamiltonian 
interfere destructively for the dark state. Using Eqs. (|73|) and (|73|) , for the 
exact expression of the dark state Eq. (|55|), one obtains 



nN Pac {0)\D$ 

^ Pac (o)x:(-i) r 



(C4) 



m=0 



/n(n — 1) • • • (n — m + 1) 



nr. 



X' 



'JV(JV - 1) • ■ • (JV - m + 1) fi n " m (c/VAO m 
JV™ (fi 2 + g 2 N) n / 2 



\n-m)\C™) 



and 



E(-i) r 



m=l 



\ 



n(n - 1) • • • (n - m+ 1) /JV(JV - 1) • • • (JV - m + 1) 



Q n - m+1 (gVN) n 
{Q 2 + g 2 N) n / 2 



(m- 1)! 

In-m)!^,^" 1 ) 



JV r ' 



»Mff)p«*(g)P?> 

re 
m=0 



'n(n — 1) • • • (n — m + 1) 



ml 



x< 



JV(JV - 1) • • • (JV - m + 1) fi n " m (c/\/JV) 



(C5) 



JV' 



(n 2 + c/ 2 jv)< 2 



|n-m)iO 



n-l 

EC- 1 ) 7 

m=0 

n n - m (gVN) 



ln(n- 1) •••(n-m) N(N — 1) ■ ■ • (JV — m) 



m! 



,m+l 



(ft 2 + ^ 2 JV)"/2 



|n-m-l)L4^), 
25 



If we set index m — > m + 1 in Eq. (|73"| ) , it is the exact same as Eq. ([73]) but 
with an opposite sign. Therefore, one obtains 

[gNa(q)p ab (q) + m Pae (0)]\^) = 0. (C6) 

Resulting in, for the exact resonant mode q — 0, the dark states \D™ =0 ) being 
the eigenstates with the null eigenvalue of the interaction Hamiltonian ([49]) . 
We note that Eqs. ([7jJ) - (|7"3|) hold exactly for the exact expression of the 
dark state Q55|) . For the approximate expression of the dark state fl56"D , Eqs. 
(HI) ~~ (HI) are & l so satisfied as long N >> n. 
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